Electromagnetic and gravitational responses and anomalies in topological insulators 

and superconductors 
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One of the defining properties of the conventional three-dimensional ("Z2-", or "spin-orbit"-) 
topological insulator is its characteristic magnetoelectric effect, as described by axion electrody- 
namics. In this paper, we discuss an analogue of such a magnetoelectric effect in the thermal (or 
gravitational) and the magnetic dipole responses in all symmetry classes which admit topologically 
non-trivial insulators or superconductors to exist in three dimensions. In particular, for topologi- 
cal superconductors (or superfluids) with time-reversal symmetry which lack 5(7(2) spin rotation 
symmetry (e.g. due to spin-orbit interactions), such as the B phase of 3 He, the thermal response is 
the only probe which can detect the non-trivial topological character through transport. We show 
that, for such topological superconductors, applying a temperature gradient produces a thermal- 
(or mass-) surface current perpendicular to the thermal gradient. Such charge, thermal, or mag- 
netic dipole responses provide a definition of topological insulators and superconductors beyond the 
single-particle picture. Moreover we find, for a significant part of the 'ten-fold' list of topological 
insulators found in previous work in the absence of interactions, that in general dimensions the ef- 
fective field theory describing the space-time responses is governed by a field theory anomaly. Since 
anomalies are known to be insensitive to whether the underlying fermions are interacting or not, this 
shows that the classification of these topological insulators is robust to adiabatic deformations by 
interparticle interactions in general dimensionality. In particular, this applies to symmetry classes 
Dili, CI, and AIII in three spatial dimensions, and to symmetry classes D and C in two spatial 
dimensions. 



I. INTRODUCTION 



The considerable recent progress in understanding 
topological insulating phases in three dimensions was ini- 
tiated by studies of single-particle Hamiltonia|HS describ- 
ing electrons with time-reversal invariance.uu In both, 
two and three dimensions, time- reversal invariant Fermi 
systems which have topological invariants of Z2 type are 
known to cxists:insulators can be classified as "ordinary" 
or "topological" by band-structure integrals similar to 
the integer-valued integrals that appear in the integer 
quantum Hall cffcct.u'Ll These invariants survive when 
disorder is added to the system. In fact, stability to dis- 
order is one of defining properties of topological insulat- 
ing phases (and also topological superconductors). The 
complete classification of topological insulators and topo- 
logical superconductors in any dimension has been ob- 
tained in Refs. || and ^, and in every dime ns i on five of the 
ten Altland-Zirnbauer symmetry classealilO of single- 
particle Hamiltonians (including some describing the Bo- 
goliubov quasiparticles of superconductors or superfluids, 
rather than ordinary electrons) contain topological insu- 
lating phases with topologically protected gapless surface 
states. 

An important question is how these various phases 
can be defined in terms of a physical response function. 
Aside from aiding in experimental detection, such def- 
initions also indicate that the phase is well-defined in 
the presence of interactions. The best studied example 
is the conventional three-dimensional ( "Z2-" , or "spin- 



orbit"-) topological insulator with no symmetries beyond 
time-reversal, which has been recently observed in var- 
ious materials including Bi^Sbi-a; alloyai3, Bi2Se3, and 
Bi2Te3oil3 Such materials support a quantized magne- 
toelectric response generated by the orbital motion of the 
electrons, i.e., the phase can be defined by the rcsppnse. 
of the bulk polarization to an applied magnetic ficld£3'E£l 
The possibility of such a bulk response was discussed 
some time ago as a condensed matter realization of "ax- 
ion electrodynamics" B3 

The first goal of this paper is to find, for all three- 
dimensional topological insulators and superconductors, 
the corresponding responses that result fcpm the coupling 
of the theory to gauge and gravitational^!] fields. The sec- 
ond goal of this paper is to understand to what extent 
the classification scheme found previously for topologi- 
cal insulators of non-interacting fermions can be stable 
to fermion interactions. This addresses the question as 
to whether certain topological insulators which describe 
distinct topological phases in the absence of fermion in- 
teractions (connected only by quantum phase transitions 
at which the bulk gap closes), can be adiabatically de- 
formed into each other when interactions are included 
(without closing the bulk gap) . We find that this cannot 
happen e.g. in symmetry classes Dili, CI and AIII in 
three spatial dimensions, and in symmetry classes D and 
C in two dimensions. More generally, in the last (more 
technical) chapter of this paper we provide an answer to 
this question in general dimensionalities for a significant 
part of the list of topological insulators (superconductors) 
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within the 'ten-fold' classification scheme, obtained for 
non-interacting particlesocjc3 In particular, we relate 
the topological features of these topological insulators to 
the appearance of a topological term in the effective field 
theory describing space-time dependent linear responses. 
Alternately, we relate these topological terms to what are 
known as 'anomalies' appearing in the theories describ- 
ing the responses. Since the 'anomalies' are known to be 
insensitive to whether the underlying fermions are inter- 
acting or not, our so-obtained description of the topo- 
logical features demonstrates the insensitivity of these 
topological insulators (superconductors) to adiabatic de- 
formations by interactions. 

The general picture emerging from the results pre- 
sented in this paper is that the topological insulators 
(superconductors) appearing in the "ten-fold list" can be 
viewed as generalizations of the d = 2 Integer Quantum 
Hall Effect to systems in different dimensions d and with 
different ( "anti- unitary" ) symmetryO properties. While 
the "teiiplbld classification scheme" was originally estab- 
lished intjij for noninteracting fermions, the characteriza- 
tion in terms of anomalies implies that this extends also 
to all those interacting systems which can be adiabati- 
cally connected to noninteracting topological insulators 
(superconductors) without closing the bulk gap. (This 
may include fairly strong interactions, albeit typically not 
expected to exceed the noninteracting bulk gap.) One 
may expect that to any of the topological insulators (su- 
perconductors) in the "ten-fold list" (viewed as general- 
izations of the Integer Quantum Hall Effect) corresponds 
a set of "fractional" topological insulators (superconduc- 
tors) not adiabatically connected to a noninteracting one, 
in analogy to the case of the two-dimensional Quantum 
Hall Effect. This includes, e.g., a recently proposed three- 
dimensional "fractional" topological insulatorn3. One ex- 
pects a description in terms of anomalies to carry over 
to all such systems and to play a role in a (future) per- 
haps comprehensive characterization of such "fractional" 
topological insulators (superconductors). In the present 
paper, however, we focus on those interacting topological 
insulators (superconductors) which can be adiabatically 
connected to a noninteracting system of fermions. 

Let us focus now on the topological insulators (super- 
conductors) in d = 3 spatial dimension (see also Table 
||) . From a conceptual point of view it is the surface re- 
sponses which are simplest to describe and they are quan- 
tized (but they may not necessarily be most easily acces- 
sible experimentally; therefore we also discuss the bulk 
responses further below): 

Charge surface response: this is, in particular, relevant 
for the ("Z2", or "spin-orbit") topological insulator which 
is time reversal invariant. Upon subjecting its surface 
to a weak time reversal symmetry breaking perturbation 
(in the zero temperature limit), the surface turns into 
a quantum Hall insulator whose electrical, surface Hall 
conductance takes on the quantized valuecil 



(a multiple of half the conductance quantum) as the 
strength of the symmetry breaking perturbation is re- 
duced to zero (always at zero temperature). Here n = 
and n =j_ .L-.for the "Z2" (or "spin-orbit") topological 
insulatoniiJE-] (in the so-called symmetry class All), in 
the topologically trivial and non-trivial phase, respec- 
tively. While the surface of Z2 topological insulators in 
class All may exhibit any odd (even) number Dirac cones 
in the topologically non-trivial (trivial) phase at the mi- 
croscopic level, only the odd-even parity, n = 1 and n = 
of that number is topologically protected. On the other 
hand, for the less familiar topological insulator in sym- 
metry class AIII, it will become evident from the results 
of the present paper that the number n in the surface 
Hall conductance (0) is not restricted to n = or 1 but 
the maximal (in magnitude) value it can take on is equal 
to the topological charge (winding number) v defined for 
the bulk topological states.E3c3 

Spin surface response: analogous effects are knownB for 
the time reversal invariant topological (spin-singlet) su- 
perconductor in symmetry class CI in d = 3 spatial di- 
mension. Subjecting its surface, as above, to a weak time 
reversal symmetry breaking perturbation (in the zero 
temperature limit), the surface turns into .what is known 
as the "spin quantum Hall insulator" E3e3. Due to spin- 
singlet pairing this superconductor has SU(2) Pauli-spin 
rotation symmetry which permits the definition of the 
'surface spin conductivity'. In particularEH, a gradient of 
magnetic field within the surface (say in the z-direction 
of spin space) leads to a spin-current perpendicular to the 
gradient (and within the surface). This defines the 'sur- 
face spin-Hall conductance' which takes on the universal 
value 



(spin) 



(2) 
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(n-times half the 'spin-conductance quantum' 
where the maximal (in magnitude) value the integer 
n can take |Qp-|is equal to the integer v from the Z- 
classificationtjEa of the topological superconductor in 
symmetry class CI) as the time reversal symmetry break- 
ing perturbation is reduced to zero. 

Thermal surface response: as we show in section [II B of 



this paper, an analogous effect occurs for the thermal re- 
sponse at the surface of the time reversal invariant topo- 
logical superconductor in symmetry class Dili in d = 3 
spatial dimensions: subjecting its surface, as above, to 
a weak time reversal symmetry breaking perturbation 
(in the low temperature limit), a temperature gradient 
within the surface leads to a heat (energy) current in the 
perpendicular direction in the surface. The so-defined 
surface thermal Hall conductance <rj y (when divided by 
temperature) tends, as will become evident from the re- 
sults of the present paper, in the zero-temperature limit 
to a quantized value 



,/(e 2 /h) 
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±c/2, where c = n/2 (3) 
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where the maximal (in magnitude) value the integer 
n can take iQp-iis equal to the integer v from the Z- 
classificationHea of the topological superconductor in 
symmetry class Dili, as the symmetry breaking pertur- 

bation is reduced to zero, [c x i7r „f ' is the thermal 



Symmetry Charge Gravitational 



Dipole 



3ft, 



conductance for a Majorana fermion when c = 1/2 (its 
central charge).] 

These surface Hall conductances provide a characteri- 
zation of the mentioned topological insulators, irrespec- 
tive of whether electron interactions are present or not: 
if we start out with a non-interacting topological insu- 
lator, one can explicitly compute the theory describing 
various space-time dependent responses. (For the ther- 
mal responses of the Dili topological supercond uctor in 
d = 3 spatial dimension, this is done in Subsec. IHB of 
this paper. For the SU(2) spin-responses of the topo- 
logical singlet superconductor in symmetry class CI this 
was done in Ref. For a significant part of the list of 
all topological insulators (superconductors) this is done 
more generally in Section ^ of this paper for all dimen- 
sionalities.) Owing to the fact that the underlying insu- 
lators are topological, the field theories for the responses 
turn out to be described by what are called anomalies. 
These anomalies describe the responses both, in the bulk 
and at the surface. The charge, spin, and thermal sur- 
face responses are examples £3 Anomalies are known to 
be insensitive to the presence or absence of interactions. 
For this reason, the (maximal value of the) dimensionless 
charge, spin, and thermal surface responses are indepen- 
dent of the strength of the interactions. These surface 
responses can only change when a bulk quantum phase 
transition is crossed (at which the bulk gap closes). 

While these surface responses are quantized and the- 
oretically useful in that they permit one to understand 
the stability of the topological insulator (superconductor) 
phases to interactions (for the cases discussed above, and 
in Section ^ for general dimensionalities) , they may not 
all be directly accessible experimentally. Therefore we 
discuss below also the various bulk responses. 

The bulk responses that we find arc of three types: 
charge response, previously shown to lead to a quantized 
E ■ B term in the ordinaxi£JZ 2 topological insulator ( "ax- 
ion electrodynamics" ) Ii3c3; gravitational response, when 
energy flows lead to an analogue of this term for gravita- 
tional-Jjields, leading to a Lense-Thirring frame-dragging 
effected when a temperature gradient is applied; and mag- 
netic dipole response, when a magnetic dipole current in- 
duced by an applied perturbation leads to an electrical 
field. A single phase may show more than one of these 
effects; for example, a phase with a conserved SU (2) spin 
current can show a non-Abclian response of this type in 
the presence of an SU{2) gauge field coupling to this 
current, but will also show a magnetic dipole response 
via its coupling to ordinary U(l) electromagnetism. We 
obtain these possible responses for each of the five sym- 
metry classes in three dimensions supporting topological 
phasesffl As in the classification in Ref. |[ the approach 
we take is based upon the surfaces of these topological 



All 

CI 

CII 

Dili 

AIII 
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TABLE I. Electromagnetic and gravitational (thermal) re- 
sponses for five out of ten Altland-Zirnbauer symmetry classes 
(All, CI, CII, Dili and AIII). The assumptions made in the 
first four classes are that (7(1) conserved currents arise from 
electrical charge and that SU (2) conserved currents arise from 
spin. *: In class AIII, the U(l) conservation law may arise 
either from charge or one component of spin. 



phases; these surfaces carry currents leading to new terms 
in the effective action of gravitational and electromag- 
netic fields. Our results for the various symmetry classes 
with topological invariants in three dimensions are sum- 
marized in Tabic Q. 

These bulk responses are "topological" to varying de- 
grees. The charge response is topological both in its spa- 
tial dependence and as a term of the effective action: 
quantization of the response is tied to quantization of 
the electrical charge and the Dirac quantization condi- 
tion. The gravitational response is topological in terms 
of the spatial dependence, but its coefficient is related 
to the mass or energy of the underlying particles and 
hence not quantized to the same degree as the charge 
response. The magnetic dipole response is not topologi- 
cal in the sense of being metric-independent, but it does 
arise from sample boundaries in the same way as the 
other responses. 

This paper is organized as follows: We begin in Sec. 
by reviewing the axion electromagnetism for the three- 
dimensional topological insulators in the spin-orbit sym- 
metry class (symmetry class All). In Sec. Ill, thermal re- 
sponse of three-dimensional time-reversal invariant topo- 
logical superconductors (such as the B-phase of 3 He) is 
discussed by exploiting a close analogy of electromag- 
netism and gravity in Newtonian approximation. In Sec. 
IV, dipole response is discussed for three-dimensional 
topological phases when at least one component of spin 
is conserved. All these responses will be discussed from 
much a broader perspective in Sec. [v]in terms of anoma- 
lies of various kind (chiral anomaly, gauge anomaly, 
gravitational anomaly) , and the descent relation relating 
these anomalies. We conclude in Sec. VI. 



II. CHARGE RESPONSES 

For an explicit example, consider a cylinder of topolog- 
ical insulator with surface Hall conductance ±e 2 /(2/i), 
defined with reference to the outward normal (see Fig. 
|l|). (Below, we choose plus sign for the surface Hall con- 



4 



ductance, by subjecting the surface to a weak external 
time-reversal symmetry source.) The motivation for con- 
sidering this example in some detail is that it will lead to 
a direct interpretation of the corresponding gravitational 
response below. The current response to an applied elec- 
trical field along the cylinder axis is 



j=j$0 where j e = ^j-E z 



(4) 



Now the magnetic field induced by this current follows 
from one of Maxwell's equations, 



(a) 



(b) 



J or j 



E or VT/T 



■CD ^S 1 



^c£- ■ 



_ 47T. 

V x B = — j, 

c 



which leads to the magnctostatic equation 



B(x) 



j(x') x 



(x-x') 

lv - v'|3 



dV. 



(5) 



(6) 



FIG. 1. Electric and thermal response of topological insula- 
tors, and thermal response of topological triplet superconduc- 
tors, in a cylindrical geometry, (a): Electric (j) or thermal 
(j T ) current driven by applied electric field (E) or thermal 
gradient (VT/T). (b): A response dual to (a) where an ap- 
plied magnetic field in z-direction induces charge polarization. 



The result for a thin cylinder is that the magnetic field 
at the cylinder axis, well away from the cylinder ends, is 
given by B = B z z with 



B, 



r(2nr)jg 47r 



2ne 2 E z 
he 



(7) 



This magnitude follows from minimizing the magnetic 
energy 



H B = 



21- 

8tt 2hc 



-E B 



(8) 



which follows from the Maxwell Lagrangian supple- 
mented with the theta term (axion term) 



Co 



9e 2 
2nhc 



-E B 



■e^ pX F^F pX (9) 



for the coupling 9 = —n. (The negative sign in this 
equation is picked out by the choice of the direction of 
the current flow around the cylinder.) 

To understand the dual response, which is an electrical 
field induced by an applied magnetic field, one needs to 
include the ends of the cylinder. Applying a magnetic 
field normal to a Hall layer increases or decreases the 
charge density depending on the direction of the field, as 
is required for the charge continuity equation to follow 
from Maxwell's equation 



<9B 
~dt 



V x E = 0. 



(10) 



Hence an applied magnetic field induces an electrical po- 
larization along the interior of the cylinder. We now turn 
to a gravitational version of the above physics, generated 
by energy flows from surface thermal Hall layers. 



III. GRAVITATIONAL RESPONSES 

A. gravitoelectromagnetism 

Our approach will be to start from the energy flow at 
surfaces of a topological phase, which is the microscopic 
source of the gravitational response. The importance of 
this response is that it is the only one that exists in the 
important symmetry class Dili, which includes super- 
fluid 3 He. We use this phase as an explicit example in 
the following. The surface Majorana mode that exists in 
this phase does not carry charge, but does carry heat, 
leading to a thermal Hall effect. Hence a temperature 
gradient applied to a cylinder leads to an energy flow 
perpendicular to the applied gradient 



Jo 



[11) 



where for future use we have treated temperature as a 
scalar potential generating a field E g = — c 2 (VT)/T with 
units of acceleration. The physical meaning of this scalar 
potential was worked out by Luttinger is his deriva- 
tion of the thermal transport coefficicntsE 2 !: in a near- 
equilibrium system, the effect of a thermal gradient is 
equivalent to that obtained from a gravitational poten- 
tial ip such that 



VV> = — , 



(12) 



where ip is the gravitational potential energy per mass, 
divided by c 2 . 

This rotational energy flow couples to the gravitational 
field at the first post-Newtonian approximation (i.e., the 
coupling is down by a factor v/c compared to the static 
gravitational effect present in the absence of the applied 
gradient). Because temperature couples to the local en- 
ergy density in the same way as an applied gravitational 
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potential, as used by Luttinger in his derivation of the 
thermal Kubo formula^ 2 ], we can view this effect sim- 
ilarly to the charge response above, as a gravitational 
"magnetic" field resulting from the energy flow that was 
induced by a gravitational "electric" field reflecting the 
temperature gradient. 

This analogy can be made precise in the ncar- 
Newtonia»n limit using the gravitoelectromagnctic 
equations^! that apply to a near-Minkowski metric. 
The relevant equation is that a mass current induces a 
gravitomagnetic field B g , defined more precisely below, 
via the equation 



V x B„ 



-47rGj r 



(13) 



Here j m is the (three-dimensional) mass current density, 
satisfying j m = j T /c 2 , and G is the effective Newton 
constant of the material. The negative sign in this equa- 
tion compared to the corresponding Maxwell's equation 
is physically significant and results from the difference 
that equal masses attract, while equal charges repel. The 
field Eg, like B g , has units of acceleration, and the grav- 
itational force on a test particle of small mass mt es t is 



F = mtest E„ 



v 

2- x B, 



(14) 



where v is the particle velocity. The factor of 2 here 
results from the spin-2 nature of the gravitational field. 

Now, by the same steps as above, there is an induced 
field along the cylinder axis 



B„ 



4ttG TalyE, 



'Jz 



(15) 



Since a^y nas the units k^T/h of a two-dimensional ther- 
mal conductivity, the ratio between B g and E g is of the 
form G(energy 2 )/(/ic 5 ), which is dimcnsionlcss (the grav- 
itational analogue of the fine structure constant that ap- 
pears in the charge case). 

The gravitomagnetic field then has exactly the same 
spatial dependence as the magnetic field in the axion case 
computed above. In particular, it is topological (e.g., the 
field at the cylinder axis does not fall off as the cylinder 
radius becomes larger) and scales with the energy flow, 
which in turn scales quadratically with the mass of the 
underlying particles. 



B. gravitational instanton term 

We now discuss the gravitational response in topo- 
logical insulators and superconductors from more for- 
mal point of view. When discussing electromagnetic re- 
sponses in topological insulators, we can couple electrons 
to an external (background) U(l) gauge field. The #-term 
in the effective action for the gauge field then results by 
integrating over the gapped electrons. In order to discuss 



gravitational and thermal responses, we can take a simi- 
lar approach: we can introduce an external gravitational 
field that couples to fermions (electrons for topological 
insulators, and fermionic Bogoliubov quasiparticles for 
topological superconductors.) By integrating over the 
gapped fermions, we obtain an effective gravitational ac- 
tion. The derivation of the effective action proceeds in a 
way quite parallel to that of the U(l) case: Indeed, both 
of them are related to a chiral anomaly as we will sec 
below. 

For topological insulators or superconductors defined 
on a lattice, it is not obvious how to couple fermions 
to gravity in a way fully invariant under general coordi- 
nate transformations. Also, there is of course no Lorentz 
symmetry on a lattice. Yet, energy and momentum are 
conserved, and one can think of introducing an external 
field which couples to these conserved quantities. The 
gravitoelectromagnctic approach discussed in the previ- 
ous subsection is based on a particular background (flat 
Minkowski metric), and is an approximation of the full 
Einstein gravity in the limit where the mass flows are 
small in some particular reference frame defined by the 
system with no thermal perturbation. 

Howevep^ all topological insulators (superconductors) 
are knowr£3 to possess a representative in the same topo- 
logical phase which is described by a Dirac hamiltonian. 
Fermions whose dynamics is described by a Dirac hamil- 
tonian can naturally be coupled to a gravitational back- 
ground field. (The theory is fully Lorentz invariant, and 
the coupling to gravity is fully invariant under general co- 
ordinate transformations, and can be described in terms 
of the spin connection.) For this reason we provide (be- 
low) a derivation of the effective action in terms of the 
Dirac representative of the topological phases. The topo- 
logical features of the effective action for the gravitational 
responses are expected to be independent of the choice of 
representative in the topological class, and thus to have a 
much more general applicability. Physically, such gravii 
tational responses describe thermal response functions.!! 2 ] 

We thus consider the following single 4x4 continuum 
Dirac model, 



H 



d 3 x ip' (—id ■ a + m{3) tp, 



(16) 



where ^ and ip represent creation/annihilation operator 
of complex fermions, and a = o\ ®cr and (3 = 03 ®uq are 
the Dirac matrices (co. 1,2.3 are standard Pauli matrices). 
(In this subsection, we use natural units, c = H = 1, and 
set the Fermi velocity to be one for simplicity.) For topo- 
logical superconductors, we need to use real (Majorana) 
fermions instead of complex fermions. 

We assume the Dirac model is in a topologically non- 
trivial phase for m > while it is in a trivial phase for 
m < 0: While this does not look apparent from the action 
in the continuum limit, when the Dirac model is derived 
from an appropriate lattice model, the sign of the mass 
does determine the nature of the phase. In the presence of 
a gravitational background, the fermionic action is given 
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byS 

S[m, i>,i>,e] = J d 4 Xy/gC, (17) 

C = $e a »vf(d» - |w M o6 E a6 ) V - rmH>, 

where /x, i>, ... = 0,1,2,3 is the space-time index, and 
a,b, . . . = 0, 1, 2, 3 is the flat index; e a M is the viclbcin, 
and u^ ab is the spin connection; E a h = [7 0) 7&] 
(See Ref. [jB] for our conventions of metric, viclbein, 
spin connection, etc.) The effective gravitational action 
W e ff [to, e] for the gravitational field is then obtained from 
the fermionic path integral 

e iVfe»[m,e] = J V e iS[m,^>,e] ^g-j 

A key observation is that the continuum Hamiltonian 
H enjoys a continuous chiral symmetry: we can flip the 
sign of mass, in a continuous fashion, by the following 
chiral rotation 

H^ = e' W V, ^ + -> tp 1 = ijj v e- i(h5/2 , (19) 

under which 

i> {id^l^ - m) i/j = (id^ - ml '((/>)) tp' , 

m'(4>) = ™e M75 = to [cos0 + ij 5 sin0] , (20) 

so that m'((f> = 0) = m and m'(0 = 7r) = —to. Since m 
can continuously be rotated into — m, one would think, 
naively, W^s [to, e] = W c ff [— to, e]. This naive expectation 
is, however, not true because of chiral anomaly. The 
chiral transformation which rotates to continuously costs 
the Jacobian J of the path integral measure, 

V[$,li\ = JVlfi,i>'}. (21) 

The chiral anomaly (the chiral Jacobian Jf) is responsi- 
ble for the 0-term. The Jacobian—/ can be computed 
explicitly by the Fujikawa methods with the result 

W e eS :=-\nJ (22) 

when to > while W® s = when m < 0. The expression 
in square brackets is the so-called Dirac. genus (see Sec- 
tion ^ below for details) which is equalEJ, by the Atiyah- 
Singer index theorem, to the index of the Dirac operator 
in the curved background. The multiplicative prcfactor 
1/2 arises because of the Majorana nature of the Bogoli- 
ubov quasiparticlcs. The index in square b:cackcts is in 
fact an even integer (by Rochlin's TheoremEj). There- 
fore, (1/2) of that expression, i.e. half the index, is an 
integer. Thus the gravitational effective action W% s in 
Eq. ( p2[) equals 9 times an integer, i.e. it is a so-called 9- 
term. Now, since 9 — > —9 under time reversal, the theta 
angle is fixed by time reversal symmetry and periodic- 
ity to either 9 — or 9 = tt. The former corresponds 



to a topologically trivial state, and 9 = ir to the topo- 
logically non-trivial state. (See, for a similar discussion 
on the derivation of the fl-term, i.e., E • B term, for the 
electromagnetic response, Ref. |26[ as well as for the non- 
Abelian SU(2) response in Ref. |29|.) Note that if instead 
we consider complex (Dirac) fermions in the background 
gravity field, the theta angle 9 is an integer multiple of 
2ir, but not of tt as in the Majorana case. 

The part of the effective action, which is not related 
to the Fujikawa Jacobian, takes the form of the Einstcin- 
Hilbcrt action W-eh = (WttG)" 1 J d A x *JgR where G is 
the effective Newton constant in the bulk of the topolog- 
ical insulator (superconductor). The gravitoelectromag- 
netism equations mentioned above can be derived from 
the effective action by taking the Newtonian limit (near 
Minkowski limit). 

To make the connection with the existence of topolog- 
ically protected surface modes we note that when there 
are boundaries (say) in the x 3 -direction at x 3 = L+ and 
at x 3 = L-, the gravitational instanton term Wg ff , at the 
non-trivial time-reversal invariant value 9 = tt of the an- 
gle 9. can be written in terms of the gravitational Chcrn- 
Simons terms at the boundaries, 

W c e s = 7 CS |**=i + - Ics\m_ , (23) 
where (i, j, k = 0, 1, 2) 

J cs = g 4^ Y\ J <Pxe iik tT(u)idjW k + ^cjiujjuj^j (24) 

with c = 1/2. This kind of relationship between the theta 
term and the Chern-Simons type term in one lower di- 
mension is a special case of the so-called descent relation 
and will be discussed further in Section This value 
of the coefficient of the gravitational Chern-Simons term 
is one-half of the canonical value (I/Att) x (c/24) with 
c = 1/2. As before, for fermions with a reality condition 
(Majorana fermions), the canonical value of the coeffi- 
cient of the gravitational Chern-Simons term corresponds 
to c = 1/2, as opposed to c = 1 for fermions without a 
reality condition (and c = n/2 for n species of fepmions 
with a reality condition.) As discussed by VolovikED and, 
Read and GreenEI in the context of the two-dimensional 
chiral p-wave superconductor, the coefficient of the grav- 
itational Chern-Simons term is directly related to the 
thermal Hall conductivity, which in our case, is carried 
by the topologically protected surface modealj. (See Eq. 
(|J) of the Introduction.) 

IV. DIPOLE RESPONSES 

A. topological singlet superconductor (class CI) 
and spin chiral topological insulator (class CII) 

The last response we consider can be measured in sys- 
tems with a conserved spin or magnetic dipolc current. 
Among the five symmetry classes which admit a topolog- 
ical phase in three-spatial dimensions, we thus focus on 
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topological singlet superconductors in symmetry class CI 
(possessing time- reversal and spin rotation invariance), 
and also on topological insulators in symmetry class CII 
(possessing time-reversal but without spin rotation in- 
variance) (see Table Q) . 

Simple lattice models of the three-dimensional topo- 
logical singlet superconductor in symmetry .class CI were 
discussed previously on the diamond latticeEEl and on the 
cubic latticecj. for which, in the presence of a boundary 
(surface), there is a stable and non- localizing Andreev 
bound state. Similar to the quantized E • B term for the 
charge response in the topological insulator, the response 
of topological singlet superconductors to a fictitious ex- 
ternal SU(2) gauge field ("spin" gauge field which cou- 
ples to conserved spin current) is described by the #-term 
at 9 =j_f in the (3+1) dimensional S77(2)' Yang-Mills 
theory.o The 0-term predicts the surface quantum Hall 
effect for spin transport (the spin quantum Hall effect), 
as already mentioned in the Introduction (Section Q). 

To detect such a quantum Hall effect for the SU{2) 
symmetric spin current requires a fictitious external spin 
gauge field, and hence one would think it cannot be de- 
tected experimentally. Nevertheless, we discuss in this 
section that the electromagnetic response carried by the 
dipole moment of the spin current can be measurable. 
(See Ref. [O for a similar discussion on the dipole re- 
sponse in a^tle-A superfluid thin film or two-dimensional 
p-wave paired states. ) 

The topological insulator in symmetry class CII (called 
"spin chiral topological insulator" in Ref. p6| ) is in many 
ways analogous to the more familiar quantum spin Hall 
effect in two spatial dimensions, but requires the chiral 
symmetry in addition to time-reversal symmetry. (For a 
lattice model of the Z2 topological insulator in symme- 



wherc a is some constant. .Here we have introduced a 
scalar field 9 ("axion" field)c3, by analogy with the lo- 
cal electromagnetic polarizability of the (All, spin-orbit) 
topological insulator, to describe the spatial location of 
the dipole current, which as before is a surface property. 
Here jf represents the a-th component of a magnetic 
dipole current of dipoles in spatial direction i. Such a 
current can generate two types of static electromagnetic 
responses: a dipole density through the continuity equa- 
tion 



daf + d t n a = 0, 
and an electrical field through the equation 



(V x E). = e ljk djE k 



JL 

4tt 



doit, 



(26) 



(27) 



where fi is the permeability of the material of inter- 
est. (One could alternately have a time- varying magnetic 
field, just as a current density can produce either a con- 
stant magnetic field or a time-varying electrical field.) 
The second response may be unfamiliar but can be de- 
rived from elementary principles; see Ref. ^2| for a dis- 
cussion of how it can be measured experimentally. Start 
from a dipole field in the lab frame; take one copy with 
the dipoles pointing along some direction n and boost 
that along v, and take another copy with the dipoles 
pointing along — n and boost that along —v. For a dipole 
density n a , this leads, in the comoving frame, to the field 
B a — (/V47r)n a , and hence 



V ■ B = ^-d a n a . 



(28) 



try class CII, See Ref. |26|). Just like an intuitive under- Tj sing tne non-relativistic Lorentz transformation law 



standing of the quantum spin Hall effect can be obtained 
by starting from two decoupled and independent quan- 
tum Hall systems with opposite chirality for each spin 
and then glue them together, this spin chiral topological 
insulator can be obtained by considering two indepen- 
dent topological insulators in symmetry class AIII. More 
general quantum spin Hall states or spin chiral topologi- 
cal insulators can then be obtained by destroying the S z 
conservation by mixing spin up and down components. 
The dipole response for class CII topological insulators, 
which we will describe below, assumes that a U(l) part of 
the SU{2) spin rotation symmetry is conserved (i.e., one- 
component of spin is conserved). However, even when 
there is no such symmetry, if mixing between two species 
is weak, we can still have such a dipole response. 



E->7(E + vxB) 



(29) 



with 7 ~ 1 leads to Eq. @, with jf = . 

Now we consider these responses for the surface spin 
current of a three-dimensional topological singlet super- 
conductor. The spin Hall current is always divergence- 
free by commutation of derivatives, 



ddt = oce ijk di (dj9d k B a ) = 0, 



(30) 



since whichever term the di acts on gives zero. However, 
the electromagnetic response can be nonzero: 

eijkdjEh = -j^daj? = ^-d a (ei mn d m 9d n B a ) . (31) 



B. magnetic dipole responses 

The spin current response at the surface of such a sys- 
tem to an applied magnetic field B via the Zceman effect 
can be written as 



ae. 



jk{dj9)dkB a , 



(25) 



There are two parts to this: one "monopole" part is only 
nonzero if d a B a ^ 0, and we therefore neglect it. There 
is also a term 



47T 



1 I mi 1 

(d a d m 9)d n B a . 



(32) 
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or below the original monopole), is 



FIG. 2. Surface of a spin chiral topological insulator (class 
CII) or topological singlet superconductor (class CI). 



C. example 

As an example, we compute this response for the case 
of a surface of a topological singlet superconductor, where 
the theta angle 9 varies as a function of the distance from 
the surface (Fig. |^). For the response to be nonzero, we 
need a = m = z, so the response is to the z component 
of the magnetic field. We get, up to a possible sign, 



(V x E) a 



~dl9d v B z 



(V x E) y = %%8a.B. 



4ir 



(33) 

For the case where 9 is first constant, then changes lin- 
early in z within a surface surface layer, and is then con- 
stant again outside this layer (Fig. |]), this response will 
occur entirely at the top and bottom surfaces of the re- 
gion of linear change. As an example relevant to possible 
experiments, we compute this response for the magnetic 
field produced by a magnetic monopole field of strength 
q m (i.e., from one end of a long magnetic dipole), sus- 
pended a distance zq above a spin Hall surface layer 
where 9 changes linearly across a thickness d. This sur- 
face layer gives two surfaces with 

(V x E) a = j? = T(3d y B z , (V x E). y = j™ = ±f3d x B z . 

(34) 

where (3 = (a/j.) / (4tv) n/d. At the top layer, the z com- 
ponent of magnetic field is, in cylindrical coordinates, 



B, 



q m zo 



(35) 



(r 2 + z 2 ) 3 / 2 ' 

which leads to a surface magnetic current of magnitude 



Je 

at the top surface. Since 
E(r) = 



3fiq m z r 

( r 2 +z 2)5/2- 



d 3 r / ( r - r/ ) X j( r 

It* — t»' 1 2 



(36) 



(37) 



we obtain that the electrical field from the top surface, 
at a height zi above the top surface (and directly above 



E,(zi) 



(9nr)dr WqmZ ° T V — 

(r 2 + 2 2 ) 5/2 r 2 + zi 2 



Evaluating this at the original height zq gives 



E z {z ) = (6ir/3q m z ) 



15z 4 



5z 3 



(38) 



(39) 



Comparing this to the case of an image charge above a 
metal, we sec that the electrical field falls off by one more 
power of height. From the above, the dipole currents are 
localized to the top and bottom surfaces of the region 
where 9 changes. The bottom surface contributes with 
an opposite sign and with z — > z + d, so we obtain 



E z (z ) 



4ir(3q n 



(z - 3 - (zo + d)~ 3 ) 



(40) 



so that for d <C z a the electric field falls off as the fourth 
power of distance. 

We can understand the scaling of the result by not- 
ing that q m divided by length cubed has units of mag- 
netic field per length; multiplying by /3 converts this to a 
2D magnetic charge current density, which has the same 
units as electric field. While the dipole response origi- 
nates in a topological phase, it is not itself "topological" 
but depends sensitively on the geometry used to probe 
it. 



V. TOPOLOGICAL FIELD THEORIES FOR 
SPACE-TIME DEPENDENT RESPONSES IN 
TOPOLOGICAL INSULATORS AND 
SUPERCONDUCTORS IN GENERAL 
DIMENSIONS FROM ANOMALIES 

The previous sections of this paper complete the list 
of the (topological) field theories describing the space- 
time dependent responses of all topological insulators 
and superconductors in three spatial dimensions (3 + 1 
space-time dimensions). In this section we will describe, 
more generally, the (topological) field theories for such 
responses in general dimensions. Most importantly, the 
main result obtained in this section is a general connec- 
tion between the appearance of such topological terms in 
the field theories for the responses and the appearance 
of what are called anomalies^ for the field theories in 
those space-time dimensions in which topological insu- 
lators (superconductors) appear. In fact, we may ask if 
the existence of a particular type of anomaly in a given 
dimension allows us to predict the existence of a topo- 
logical insulator (superconductor) of the 'ten-fold' clas- 
sification in that dimension. The answer to this ques- 
tion is affirmative. As we demonstrate below, a large 
part of the 'ten-fold' classification can be derived from 
the existence of the known anomalies in corresponding 
quantum field theories in space-time. This can then be 
thought of as yet another derivation of the 'ten- fold' clas- 
sification, in addition to the previously known derivations 
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such as that based on Anderson localization at the sam- 
ple boundaries!!, and K-theocyJ (as well as a later point 
of view based on D-branesMr 7 !). Moreover, and most im- 
portantly, the appearance of an anomaly is a statement 
about the respective quantum field theory (of space-time 
linear responses) independent of the assumption of the 
absence of inter-particle interactions. Thus, anomalies 
provide a description of topological insulators (supercon- 
ductors) in the context of interacting systems. 



A. Topological insulators (superconductors) in the 
two complex symmetry classes A and AIII from 
anomalies in the gauge field action 

1. the integer quantum Hall effect (class A) 

Let us begin by describing the topological field theo- 
ries describing the space-time dependent responses of the 
two "complex" symmetry classes, classes A as^Ani in 
the Cartan (Altland-Zirnbauer) classification.&ll§E3 This 
includes the most familiar example, namely the integer 
quantum Hall insulator (IQH), belonging to symmetry 
class A. In both symmetry classes, A and AIII, there 
has to exist a conserved U(l) charge (particle number). 
This is the electromagnetic charge, since these symme- 
try classes can be realized as normal electronic systems 
(as opposed to superconducting quasiparticle systems)c3. 
Therefore we can minimally couple these topological in- 
sulators to an external U(l) gauge field. The field the- 
ory describing the space-time dependent linear responses 
of the topological insulator can then be obtained by in- 
tegrating out the gapped fermions. The fact that the 
underlying insulator is topological is reflected in the fact 
that the effective action for the external U{1) gauge field, 
describing the electromagnetic linear responses, contains 
a term of 'topological origin', such as e.g. a Chern- 
Simons- or a 0-term, or corresponding higher dimensional 
analogues of these terms (see below for more details). 

In turn, the presence of terms of topological origin 
in the so-obtained effective action for the external U(l) 
gauge field are closely related to the presence of a so- 
called anomaly. To sec how an anomaly for the theory 
of the external U(l) gauge field can actually predict the 
presence of a topological phase, let us consider first, as 
the simplest example, the IQH insulator in d = 2 spatial 
dimensions - symmetry class A. (The space-time dimen- 
sion is thus D = 2 + 1). In fact, let us first focus attention 
on the theory of the sample boundary (the edge state), 
which has d = 1 spatial dimensions. It is known (see be- 
low) that the effective theory for the linear responses of 
the U(l) gauge field in D = 1 + 1 space-time dimensions 
(i.e. of the edge state) can have what is called a "gs 
anomaly" since the space time dimension D is evenJ 
The presence of this anomaly simply means that U(l) 
charge conservation is spoiled by quantum mechanics. In 
the condensed matter setting of the IQH insulator the 
meaning of this anomaly is that the system (i.e. the 



edge) in D = 1 + 1 space-time dimensions, exhibiting 
the anomaly, does not exist in isolation, but is necessar- 
ily realized as the boundary of a topological insulator in 
one dimension higher. In this case, the breakdown of 
the conservation law of U{1) charge conservation at the 
boundary simply means that the current "leaks" into the 
bulk. Thus, in the condensed matter setting, the pres- 
ence of the anomaly in the theory at the boundary is not 
something abnormal, but it is a physical effect: it is the 
integer quantum Hall effect. As we will discus shortly 
below, the same reasoning applies to all even space-time 
dimension, D = 2k. Consequently, we see that the pres- 
ence of a U(l) gauge anomaly predicts the presence of a 
topological insulator in one dimension higher. I.e., this 
predicts the presence of a topological insulator in sym- 
metry class A in D = 2k + 1 space-time dimensions, in 
agreement with the 'ten-fold' way classification. 



three-dimensional insulator (superconductor) in 
symmetry class AIII 



Let us now consider the topological insulator (super- 
conductor) in the other complex symmetry class, class 
AIII, in d = 3 spatial dimensions. Again, the space- 
dimension D = 3 + l = 4is even. It is known (see below) 
that in all even space-time dimensions the effective action 
for the space-time dependent U(l) gauge field may also 
possess a different anomaly [in contrast to the discussion 
in subsection above], often referred to as the "chimj (or: 
axial) anomaly in a background U(l) gauge field" £j The 
meaning of such an anomaly can be explained using Eq. 
( p6| ) below: the so-called axial (or: chiral) U(l) current 
J^{x) is not conserved in the presence of a background 
U(l) gauge field, i.e. D^J^{x) ^ 0, where denotes 
the covariant derivative in the presence of a background 
gauge field. In the simplest case of a single copy of a mas- 
sive Dirac fermion (mass m), this covariant derivative of 
the current is given by Eq. (^6|) below. As displayed in 
this equation, there are two sources of the lack of con- 
servation: (i) a finite mass m ^ and (ii) the extra 
"anomaly" term A2n+2 (to be discussed in more detail 
below), which represents the breaking of the conserva- 
tion of by quantum effects.E3 Now, as discussed in 
Ref. |2(| the presence of a "chiral (or: axial) anomaly in 
a background U(l) gauge field" implies directly the pos- 
sibility of having a non- vanishing #-term when deriving 
the effective action for the external U{1) gauge ficldEil. 
(The (9-angle is fixedEj to 6 — n by a discrete symmetry, 
which is the chiral symmetry for symmetry class AIII.) 
Thus, the presence of a "chiral (or: axial) anomaly in a 
background U(l) gauge field" in D = 2k space-time di- 
mensions signals the existence of a topological insulator 
in this space-time dimension through the appearance of 
a (9-term in the (topological) field theory for the linear 
responses. 
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3. anomaly polynomials and descent relation 

Observe that above we have used anomalies of two 
kinds, and we used them in two different ways: 

(i) : In case 1 there was an anomaly in the theory of the 
responses at the boundary [which had D = (d — 1) + 1 
space-time dimensions]. In this case the anomalous the- 
ory (i.e. the one at the boundary) was gapless (critical); 
we refer to this situation as a gauge anomaly (i.e., non- 
conservation of the U(l) charge in question). The pres- 
ence of this anomaly implied the existence of a topologi- 
cal insulator in one dimension higher, i.e. in D' = d + 1 
space-time dimensions. The responses of this topologi- 
cal insulator are described by an effective Chern-Simons 
action for the U(l) gauge field in D' = d + 1 space-time 
dimensions. [See also Eq. ji^)]. 

(ii) : In case 2 there existed an anomaly in the massive 
bulk theory in D = d + 1 space-time dimensions. This 
was a chiral anomaly (referring to the violation of the 
conservation of the global axial U(l) current J§) in the 
background of a non- vanishing U ( 1 ) background gauge 
field. 

There are important relationships between the follow- 
ing different anomalies, 

• the U(l) gauge anomaly in D = In, 

• the Chern-Simons term (i.e., parity anomaly) in 
D = 2n + 1 and 

• chiral anomaly in the presence of a background 
gauge field in D = 2n + 2, 

which can be summarized, in terms of the_so-called de- 
scent relation of the "anomaly polynomial" eft Let us now 
explain this relation. 

As mentioned above, it is known that in even space- 
time dimensions D = 2n, there is a U(l) gauge anomaly. 
If there is a gauge anomaly, the (Euclidean) effective ac- 
tion lni?L4] in the presence of the gauge field A is not 
invariant under a gauge transformation A — > A+v. Thus 
we can write 



S v In Z [A] = 2m I ^{v,A,F). 



(41) 



M 2 „ 



where the variation S v is the gauge transformation in 

question, and fi^n i s a 2n-form built from the con- 
nection 1-form, A = A^dx^, its field-strength 2- form, 
T = {\/2)F^dx^dx u , and the variation v = v^dx^ of 

the gauge field. (By definition, i s linear in v. The 
integral is taken over the physical D = 2n dimensional 
(Euclidean) space-time M2„.) Now, the descent rela- 
tion tells us that fijn can t> e derived from the so-called 
anomaly polynomial £l2n+2{J r ), which is a 2n + 2-form 
built from the curvature 2-form J- , with the aid of yet 
another 2n + 1-form f^n+ii by 



O 



2n+2 



dn 



(0) 

2n+l> 



^v^"2n J \-\ — d^l 



(1) 
2n ■ 



(42) 



I.e., fl2n+2 is closed, and gauge invariant, and hence can 

be written as a polynomial in J 7 . Here &2n+i (-4) ~F) 1S i* s 
corresponding Chern-Simons form. 

There is a simple closed form expression for the anomaly 
polynomial fi2n+2 which is given by 



tt D {F) =ch(JF)| B 



(43) 



Let us explain the notation: ch(J r ) is the following power 
series (" characteristic class" ) constructed from the field- 
strength two-form T and is given by 



Ch( ^= r+ 2^-2W 



tr T l 



(44) 



This expression is written for the general case of a gauge 
field transforming in a r-dimcnsional irreducible repre- 
sentation of a (possibly non-Abelian) gauge group, where 
tr denotes the trace in this representation. Observe that 
ch(J r ) consists of a sum of different p- forms with different 
p where p = even. The notation • • • |_d in Eq. ( [i"3f ) means 
we extract a D-form from ch(J r ). 

While up to this point the differential forms f4n+i 
and fl2n+2 appear to have been introduced solely to ex- 
press the D = 2?7-dimcnsional gauge anomaly in terms of 
other objects, they themselves are known to be related to 
other types of anomalies: the Chern-Simons form fi^n+i 
represents an anomaly in a discrete symmetry (parity or 
charge-conjugation symmetry, depending on dimension- 
ality) discussed in moife-detail in subsection V A4| , below, 
and fl2n+2 representso the chiral anomaly in the pres- 
ence of a background gauge field, discussed in subsection 
V A 2 above. The integral of ^2n+2 over D = 2n + 2 di- 
mensional space-time, on t he oth er hand, represents the 
9 term (see also subsection VA5 below). 



4- the Chern-Simons term 



The integral of n ( ° ] l+1 (A, J 7 ) over D 



2n 



1- 



dimensional space-time is the Chcrn-Simon type action 
for the gauge field A, and represents, as already men- 
tioned, an anomaly in a discrete symmetry: the parity or 
charge-conjugation anomaly. 

In turn, the presence of such a Chern-Simons term 
in the effective (bulk) action for the gauge field A in 
D = 2n + 1-dimensional space-time signals the presence 
of a topological phase: when there is a boundary in the 
system, the integral of the Chern-Simons term is not in- 
variant on its own; rather, upon making use of the de- 
scent relation Eq. (E3), one obtains 



(0) 

2n+l 



dn 



ID 



Mo 



OMo 



^•(45) 



This is something we are familiar with from the physics 
of the quantum Hall effect: the presence of the boundary 
term J gM2 i appearing on the right hand side of 
Eq. (^5|) signals the presence of edge mode. In turn, 
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as we have seen in subsection V A 2 , the gauge anomaly 
in D = 2n dimensional space-time which is represented 
by the integral over O^n ; itself signals the presence of a 
topological phase in D = 2n + 1 space-time dimensions, 
i.e. in one dimension higher. 



B. Topological insulators (superconductors) in the 
remaining eight 'real' symmetry classes from 
gravitational and mixed anomalies 



1. Gravitational anomaly and axial anomaly in the 
presence of background gravity 



the theta-term 



The integral of the anomaly polynomial ^2n+2 over 
D = 2n + 2 dimensional space-time is the #-term and 
represents a chiral anomaly in the presen ce of a back- 
ground gauge field (discussed in subsection V A2 above). 
Again, to be more explicit, in the presence of such an 
axial anomaly, the axial current J§{x) (which in the 
present case is an axial U(l) current) is not conserved: 
D^J^(x) ^ where is the covariant derivative in the 
presence of the gauge field. For a single copy of a massive 
Dirac fermion, it is given by 



Dfj,J^(x) = 2im?/>72r l +i^ + 2iA 2n +2(x), (46) 

where the first term represents the explicit breaking of 
the chiral symmetry by the mass term, whereas the 2nd 
term represents the breaking of the chiral symmetry by 
quantum effects. A 2n +2 quantifying the breaking of the 
axial current conservation by an anomaly is essentially 
identical to VL2 n +2, and given by removing all dx^ which 
appear in the differential form VL2 n +2- 

Just as it was the case for the Chern-Simons term, the 
presence of such a 0-term in the effective action for the 
gauge field signals the presence of a topological phase. In 
particular, the descent relation tells us that 



2n+2 



M 2 



dfl 



M 2 



(0) 

2ri+l 







(0) 

2ri+l 



(47) 



This is, again, something we are familiar with from the 
physics of the three-dimensional topological insulator in 
class AIII, which is described by the 0-term (the axion- 
term). In the presence of a boundary dM2 n +2i such a 
topological state supports boundary degrees of freedom, 
as signaled by the boundary term J dMi} i ^2n+i which 
is a Chern-Simons tcrm.El 

Let us summarize: in order to derive the existence of 
topological phases in symmetry class A and AIII, we start 
from the anomaly polynomial f22n+2- Then the terms 
J M 2 ^2n+2 and J M i f4n+i are the effective actions 
for the (topological) field theory of the space-time linear 
responses for the gauge field for the topological phases in 
class AIII (D — 2n + 2) and A (D = 2n+l), respectively. 



For the remaining eight "real" of the ten symmetry 
classes, having a conserved U(l) quantity is less trivial. 
Classes AI, All, and CII are naturally realized as a nor- 
mal (as opposed to superconducting) electronic system, 
and thus for these there is a natural notion of a conserved 
U(l) quantity (the electrical charge). One realization of 
the BDI symmetry class, which is only parto of the en- 
tire symmetry class, can also be considered to have a 
conserved U(l) quantity and we consider this realization 
in this subsection. On the other hand, classes D, Dili, C 
and CI are naturally realized as BdG systems. While for 
classes C and CI, SU(2) spin is conserved (so a conserved 
U(l) charge exists), for classes D and Dili, there is no 
conserved U(l) quantity at all. 

Since for the latter four of eight real symmetry classes 
(D, Dili, C, CI) we cannot rely on a conserved U(l) quan- 
tity to describe these topological phases, it is not possi- 
ble to couple these systems minimally to a U(l) gauge 
field. However, it is natural to consider a coupling of 
these topological phases to gravity. Let us focus first on 
topological insulators (superconductors) with a integer 
topological charge, Z, but not on those with a binary 
topological charge, Z2. For now we also do not consider 
topological insulators/superconductors with a 2Z charge. 

An analogue of the U(l) gauge anomaly, which we have 
described in section V A 2 at the boundary (of space-time 
dimension D = 2n) of topological phases in symmetry 
class A is the gravitational anomaly. It corresponds to 
the breakdown of energy-momentum conservation, and 
when it happens, it must be realized in a system which 
represents the boundary of a topological phase in one di- 
mension higher (in ana logy to the case of a C/(l) gauge 
anomaly, section V A2| ) . We refer to this anomaly also as 
a 'purely gravitational anomaly'. In the following we will 
show that one can predict the appearance of the topolog- 
ical phases in symmetry classes D, C, Dili, CI [i.e. those 
without conserved U(l) charge] from the presence of a 
purely gravitational anomaly which appears in the field 
theory for the gravitational (or: thermalEj) responses. 

Finally, we will need to discuss the still remaining sym- 
metry classes AI, BDI, All, and CII. Topological insula- 
tors (superconductors) in these symmetry classes can be 
coupled to both, a U(l) gauge fieldo as well as a gravita- 
tional background. We will show that the field theories 
for the space-time dependent linear responses for these 
topological insulators possess a so-called mixed anomaly. 
Indeed, we will show that the appearance of a mixed grav- 
itational and electromagnetic axial anomaly signals the 
existence of topological phases in these symmetry classes. 
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2. Topological insulators (superconductors) in symmetry 
classes D, C, Dili, CI, from the purely gravitational anomaly 

As mentioned earlier in this paper, each topological 
insulator (in aay dimension) has a Dirac Hamiltonian 
representative. E3 We can consider the coupling of this 
Dirac theory to a space-time dependent gravitational 
background. Upon integrating out the massive fermions, 
we obtain an effective gravitational action in D space- 
time dimensions. If there is a gravitational anomaly, 
the (Euclidean) effective action \nZ[e,ui] in the presence 
of the gravitational background is not invariant under a 
general coordinate transformation x^ — > x^ + e M , where 
e is the vielbein and ui is the spin-connection one-form. 
I.e., 



S v In Z[e, ui] = 2iri 



(48) 



Mr 



where S v represents an infinitesimal SO(D) rotation, un- 
der which ui, the spin-connection 1-form ui, is transformed 
as w — > ui + v; (v, u>, 1Z) is a D-form related to 
the gravitational anomaly. In complete analogy to the 
case of the gauge anomaly discussed above, Q^\v,u,TZ) 
can be derived from a corresponding anomaly polyno- 
mial VL D+2 (TZ) (see Eqs. ([54j), ( |55| ) below) through its 

Chern-Simons form n^ +1 (ui,TZ), by using a descent re- 
lation which takes a form identical to Eq. (^). Thus, 
once the existence of the (purely) gravitational anomaly 
is known for a given dimension D, it predicts the pres- 
ence of topological phases in D + 1 and D + 2 dimensions, 
using the same logic as in the gauge field case above. 

Now, according to Ref. |55], a purely gravitational 
anomaly can exist in 



D = 4k + 2 (d = 4k + l). 



(49) 



Thus, breakdown of energy-momentum conservation due 
to quantum effects can occur in these dimensions. As in 
the case of symmetry class A, discussed above, we take 
this as evidence for the existence of a topological bulk in 
one dimension higher, i.e. in space-time dimensions 

D = 4k + 3 (d = 4k + 2). (50) 

This thus predicts the appearance of topological phases 



class D [d = 2), class C (d = 



(51) 



as well as all the other higher-dimensional topological 
phases that we can obtain from these by Bott periodicity. 
(These are colored red in Table 0.) 

On the other hand, there is an analog of the "axial 
anomaly in the presence of a backgr ound gauge field" 
which we discussed in section V A2 in the context of 
symmetry class AIII in D = 2n space-time dimensions. 
This analog is the "axial anomaly in the presence of a 
background gravitational field". If only a background 



gravitational field is present, this anomaly exists in space- 
time dimensions 



1). 



(52) 



D = 4k (d = 4k 

This covers symmetry classes 

class Dili (d = 3), class CI (d = 7), (53) 

as well as all higher-dimensional topological phases that 
we can obtain from these by Bott periodicity. (These are 
colored blue in Table ||.) 

The anomaly polynomial related to the gravitational 
anomalies is known explicitly. It can be written as 

n D=ik =A{n)\ D (54) 

where A(1Z) is the so-called Dirac genus given by@ 



A(K) = 1 + 



1 1 



(4tt) 2 12 



trK z 



(4tt) 2 



(trTZ 2 ) 2 + — trft 4 

288 V 7 360 



Here 1Z is the D x D matrix of two-forms 



(55) 



(56) 



where Rap^ u is the usual Riemann curvature tensor, and 
the trace refers to the D x D matrix structure. This 
defines, by the descent relation [which takes a form iden- 
tical to Eq. ([42])], the differential forms ^ < ^_ 1 and 
As before, the notation A(1Z)\d extracts a D-form from 
A(JVj. It is obvious from (]5q) , that the anomaly polyno- 
mial exists only for D — 4k because Eq. (fjq) is a func- 
tion of 1Z 2 . (Note that the descent relation Eq. ( [42| ) then 
implies the existence of a purely gravitational anomaly 
Q^ +2 (1Z) in D = 4k + 2 space-time dimensions, in agree- 
ment with Ref. B5l) 



3. Topological insulators (superconductors) in symmetry 
classes AI, BDI, All, C II from the mixed anomaly 

Before proceeding let us briefly summarize the previ- 
ous subsection: by considering various anomalies related 
to gravity, we can predict the integer topological phases 
in the BdG symmetry classes D, Dili, C, and CI. (As 
mentioned above, for a moment, we do not consider topo- 
logical phases with Z2 or 2Z topological charges). On the 
other hand, we have so far not covered the description of 
topological insulators in symmetry classes AI, BDI, All, 
and CII in terms of anomalies. 

So far, we have considered for the 'real' symmetry 
classes only those anomalies which involve solely gravity. 
Since the (gapped) topological insulators in symmetry 
classes AI, BDI, All, and CII, also possess a conserved 
U{1) chargeE-3, we can couple those to both, a £7(1) gauge 
field as well as a gravitational background. Therefore, it 
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where 0,^ (v, A, u), T, TZ) can be deriv^imm an associ- 

reads,BB 



Sl D (TZ,F) = (ch(T)A(TZ)) \ D 



(58) 



As the right hand side is given simply given by the prod- 
uct of the anomaly polynomials for a gauge field [Eq. 
(p4|)] and gravity [Eq. (p5[)], by switching off either TZ or 
JF , we recover the results discussed in the previous subsec- 
tions: for all even spacetime dimensions D = d + 1 = 2k 
(k = 1,2,...) we obtain a non- vanishing anomaly poly- 
nomial £Id(1Z = 0,J-) = r2n(J-"), which we have already 
used to predict topological insulators/superconductors in 
class A (D = 2k + 1) and AIII (D = 2k). For space-time 
dimensions D = d+1 = 4k (k = 1, 2, . . .) we obtain a non- 
vanishing anomaly polynomial £lrj(TZ,J- = 0) = Hd(TZ), 
which we have already used to predict topological insula- 
tors/superconductors in class Dili (D = 4 + 8fc) and CI 
(D = 8 + 8k). 

On the other hand, while the anomaly polynomial 
Hd(TZ, J- = 0) = SId(JZ) vanishes in D = 4fc + 2 dimen- 
sions, the one obtained from Eq. (|58|), namely Qd(7Z,T), 
is non-vanishing in these dimensions. 



TABLE II. Topological insulators (superconductors) with 
an integer (Z) classification, (A): in the complex symmetry 
classes, predicted from the chiral U(l) anomaly, and (B): in 
the real symmetry classes, predicted from the gravitational 
anomaly (red, • • ■^), the chiral anomaly in the presence of 
background gravity (blue, ■•• ), and the chiral anomaly in 
the presence of both background gravity and (7(1) gauge field 
(green, • • •*). 



is natural to consider an anomaly which occurs in the 
presence of both, a background gauge and a background 
gravitational field. 

As it turns out, even in the presence of both gauge 
and gravitational fields, the structure of the anomaly is 
similar to the one discussed so far: the non-invariance 
of the effective action under a gauge transformation or 
coordinate transformation can be expressed as 

6 v \nZ[A,e,w] =2ni [ a < ^ > (v,A,u,J 7 ,n), (57) 

e Q,p(v,A,cj, T, TZ) can be 
ated anomaly polynomial which 



As before, the anomaly polynomial itself is related to 
a 'chiral anomaly in the presence of both gauge field and 
gravity' of the massive bulk system in D = 4k + 2 space- 
time dimensions, D^J^(x) = 2imtpjD-i'>P + 2iAd{x), 
where Ad(x) is given in terms of SI ^ (T 7 ^, . For this 
reason, one predicts an additional topological insulator 
(superconductor) in t hese s pace-time dimensions (be- 
sides the one of Sec. V A 2 ). Therefore, one predicts 
the occurrence of topological phases in spatial dimensions 
d = 9, (d = 1) and d = 5, 

class BDI (d = 9, (d = 1)), class CII (d = 5), (59) 

as well as of all higher-dimensional topological phases 
that we can obtain from these by Bott periodicity.EEl 
(These are colored green in Table [0]). Indeed, for class 
BDI and CII, we can realize these symmetry classes as 
a normal (i.e. not supercondu&ting) system, and hence 
they have a natural U(l) chargeE-1 The effective topolog- 
ical field theory for the space-time dependent linear [elec- 
trical and gravitational (thermal)] responses possesses a 
term of topological origin of the form J fijj (TZ, J-) , where 
D = 4k + 2. 

Moreover, it turns out that a descent relation which is 
identical in form to Eq. ([l2|) also holds for the 'mixed' 
anomaly polynomial defined in Eq. ([58|). Therefore, the 
space-time integral of the Chern-Simons form n^ +1 of 
^4fc+2, which is obtained from fl4k+2 by using the descent 



relation, dQ 



(0) 
4fe+l 



describes the term of topologi- 



cal origin in the effective action for the linear responses in 
D = 4k + 1 space-time dimensions. This corresponds to 
a "mixed anomaly" fij^ in the corresponding boundary 
theory in 4k space-time dimensions. For this reason, one 
predicts the occurrence of additional topological insula- 
tors in spatial dimensionalities d = and d = 4 (besides 
the ones in Sec. V A 1), for the two symmetry classes 



class AI (d = 0), class All (d = 4), 



(60) 



as well as for all their higher dimensional equivalents ob- 
tained from the Bott periodicity (These are colored ma- 
genta in Table ||). 



4- Atiy ah- Singer Index Theorem 

For all the symmetry classes with chiral symmetry the 
hamiltonian can be brought into block off-diagonal forma. 
Above we have discussed all symmetry classes of this form 
which possess topological insulators with a Z classifica- 
tion (i.e. AIII in D = 2n, Dili in D = 4 + 8k, CI in 
D = 8 + 8k, CII in D = 6 + 8&, BDI in D = 10 + 8k). 
A Dirac hamiltonian H. with chiral symmetry possesses 
an index, and the Atiyah-Singer Index TheoremL-3 relates 
the integral of the anomaly polynomial discussed above 
to this index through the formula 



index (H) 



n D {TZ,T) 



(61) 



Mr 



14 



where (1Z, J-) is the most general anomaly polynomial, 
as defined in Eq. (^38j) above. Here, the Dirac hamiltonian 
H refers to the hamiltonian in a gravitational background 
and a background (Abelian or non-Abelian) gauge field. 
The index index("H) is by definition an integer. We note 
that it is because of this theorem that the space time 
integral of the anomaly polynomial represents a 9 term 
for the theory of the space time dependent linear gauge 
and gravitational responses, and that the 9 terms only 
occur for symmetry classes possessing a chiral symmetry. 



5. Global gravitational anomalies 

The discussion that we have given so far for the con- 
nection between anomalies and topological insulators and 
superconductors in "the primary series" (those located in 
the diagonal of the periodic table and characterized by an 
integer topological invariant) can be extended to some of 
the "first and second descendants" (the topological insu- 
lators and superconductors in the same symmetry class, 
but in one and two dimensions less than the one with a 
Z invariant; these are each characterized by a Z2 invari- 
ant). We propose that for these we need to use so-called 
global anomalies, instead of the so-called perturbativc 
anomalies that we have made use of in this section. Such 
anomalies do not affect infinitesimal, but large (of order 
one) symmetry transformations. 

It was found in Ref. |55| that global gravitational 
anomalies can exist, given certain assumptions are sat- 
isfied, (i): in D = 8k, in (ii): D = 8k + 1 and (iii): in 
D = 4fc + 2 space time dimensions. If so, then follow- 
ing the same reasoning as above, the presence of these 
anomalies would indicate the existence of a topological 
insulator in one dimension higher (of which the anoma- 
lous system is the boundary). This would then indicate 
the existence of topological insulators (superconductors) 
in space time dimensions (i): D = 8fc + l, (ii): D = 8k + 2 
and (iii): D = 4fc + 3 [corresponding to spatial dimen- 
sions (i): d = 8k, (ii): d = 8k + 1 and (iii): d = 4k + 2]. 
Indeed, there exist Z2 topological insulators in these di- 
mensions (Table Q). More precisely, there exist two Z2 
topological insulators in these dimensions and at this 
point we have not yet explored in detail which of the 
two (or if both) could be related to this global gravita- 
tional anomaly. Moreover we note that there also exist 
other (i.e. not gravitational) global anomalies, and we 
propose that the other, so far not yet covered Z2 topo- 
logical insulators can be obtained from considering these 
other global anomalies. 

We end by mentioning that the notions presented in 
this Section (Section |v|) may also be further supported 
by the-cpsmcction with the tenfold classification of D- 
branesEJO: In the D-brane realizations of topological 
insulators and superconductors, massive fermion spectra 
arise as open string excitations connecting two D-branes, 
which are in one-to-one correspondence with the Dirac 
representative of the ten-fold classification of topological 



insulators and superconductors, and come quite naturally 
with gauge interactions. The Wess-Zumino term of the 
D-branes gives rise to a gauge field theory of topological 
nature, such as ones with the Chern-Simons term or the 
#-term in various dimensions. 



VI. CONCLUSIONS 

There are various important future research directions 
in the field of topological insulators and superconduc- 
tors. Let us mention two here. One is the search for 
experimental realizations of the topological singlet and 
triplet superconductors in three spatial dimensions, be- 
sides the B phase of the 3 Hc supcrfiuid. Given how fast 
experimental realizations of the quantum spin Hall ef- 
fect in two spatial dimensions and the Z2 topological in- 
sulators in three dimensions have been found, one may 
perhaps anticipate a similar development for these three- 
dimensional topological superconducting phases. No- 
tably, Cu x Bi2Sc3, which arises from the familiar three- 
dimensional topological insulators Bi2Se3, was found to 
be superconducting at 3.8 K.E3 Subsequent theoretical 
work proposed that this superconducting phase should 
be a topological superconductor E3 The various linear re- 
sponses discussed in this paper, as summarized in Table 
[jj may become helpful in the search for, and for identify- 
ing various such topological phases. 

Another important issue is to complete the study of 
the effect of interactions for the symmetry classes so far 
not yet included in the discussion given in Section 0. 
(These include, in general dimensionalities the topologi- 
cal insulators (superconductors) with a 2Z classification, 
as well as the majority of those with a Z2 classification.) 
Moreover, this includes the case of symmetry class BDI 
in d = 1 spatial dimension (recall also Rcfs. [5^ and [56| ), 
discussed in the work of Refs. ^9]-^. Further impor- 
tant outstanding questions concern possible topological 
phases (besides superconductors) which may arise from 
interactions rather than from band effects. How can 
one describe "fractional" versions of the topological in- 
sulators (superconductors) ,c3 and how can_pnc classify 
bosonic systems such as e.g. spin systems ?E3 Clearly, in 
order to address any of these interaction-dominated is- 
sues one cannot rely on a topological invariant defined 
in terms of single-particle Bloch wavefunctions. Rather, 
a definition of topological quantum states of matter in 
terms of responses to physical probes is necessary. In 
this paper we have developed a description of this type 
for all topological insulators in three spatial dimensions, 
and for a significant part of the topological insulators in 
general dimensions. From a conceptual point of view the 
gravitational responses are the most fundamental ones, 
in that they apply to alLtopological insulators. Owing 
to Luttinger's dcrivatior£3 of the thermal Kubo formula, 
these correspond physically to thermal responses func- 
tions. 
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we will call below a "chiral (7(1) gauge anomaly in the 
presence of a background gravitational field" . 
The appearance of this term for the non-Abelian gauge 
group SU(2) was first pointed out in the context of topo- 
logical insulators (superconductors) in Ref. for the spin- 
singlet topological superconductor in symmetry class CI in 
d — 3 spatial dimensions. 

For class BDI there exist two distinct physical realiza- 
tions, one as ("spinless" time-reversal invariant) supercon- 
ductors, and one as normal (non-superconducting) elec- 
tronic systems. Without considering interactions there is 
basically no difference between the two, except that the 
number of species of Majorana fermions is even in the lat- 
ter case, where a pair of Majorana fermions is thought 
to be combined into a complex fermion, carrying a U(l) 
charge, or particle number. The discussion of anomalies, 
considered in the current Section of this article, is aimed 
at the discussion of interacting theories (as explained, e.g., 
in the introduction). Now, when inclusion of interactions is 
considered, the two above-mentioned realizations of sym- 
metry class BDI behave very differently. Obviously, in the 
latter (normal, non-superconducting) realization, the in- 
teractions are to respect the U(l) symmetry, whereas in 
the former (superconducting) realization, there is no such 
constraint on the form of the interactions. In this Subsec- 
tion we will consider solely the latter realization. In this 
case, there is thus always a conserved U(l) quantity. The 
former (superconducting) case was discussed recently in 
Refs. p9j-pl|. (For similar methods applied to gapped spin 
chain, see, e.g., Ref. ^.) At present, we do not have an 
understanding of that case in terms of anomalies. We hope 
to be able to address this case in future work. 
For symmetry class BDI recall the comment in Ref. |f^. 
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For class BDI in d = 1 spatial dimensions, the mixed 
anomaly polynomial in the corresponding space-time di- 
mensionality D — 2 is simply equal to the anomaly poly- 
nomial for the U(l) gauge anomaly, discussed above (de- 
scribing the field strength in D = 2). The lowest spatial 
dimension d = 1 behaves thus differently from all other 
dimensions d — 8k + 1 (k > 1) related by Bott periodicity, 
in which an independent mixed anomaly polynomial exists. 
For this reason we have denoted d = 1 in parenthesis. 
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